A Solution Technique for Quantum Mechanical 
Differential Equations Using Multiple Complex 

Planes. 



Robert Ducharme 
December 11, 2012 



Th 1 151 Fairhills Dr., Ypsilanti, MI 48197 

O f E-mail: ducharme01@comcast.net 

Abstract 

It is shown fields that cannot be represented over one complex plane 
can be further decomposed for representation over multiple complex planes. 
This finding is demonstrated here by solving of the Schrodinger equation 
for the hydrogen atom in a complex space containing three complex planes. 
The complex coordinate system is generated from real coordinates using 
an isometric transformation. One plane is applied to mix energy and time; 
the other two planes are used to represent the z-component of angular mo- 
mentum of the electron. The eigensolutions of the Schrodinger equation 
are shown to be holomorphic in the complex planes. 

1 Introduction 

Most practical applications of conformal mapping [T] have been limited to one 
complex plane. Liouville's theorem [5] in fact shows that higher dimensional 
conformal maps are possible but must be composed of translations, similarities, 
orthogonal transformations and inversions. The purpose of this paper is to show 
problems that cannot be represented over one complex plane can generally be 
decomposed for representation over multiple complex planes. Candidates for 
decomposition in to multiple complex planes may include problems that have 
more than two dimensions in real space or two-dimensional problems that fail 
to satisfy the Cauchy-Riemann equations. The conformal transformations to 
be presented in this paper are isometric and therefore fall under Liouville's 
theorem. 

An hydrogen atom consists of a single electron with a spatial displacement Xi 
(i = 1,2,3) from the center of mass of the atom. The source of the Coulomb field 
binding the electron to the atom is the much heavier proton. Both the electron 
and the proton may be assumed to exist at the same time t. The objective of 
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section 2 of this paper, taking the Schrodinger equation of the hydrogen atom 
as an illustrative example [3], is to introduce an isometric conformal mapping 
from a real (xi,t) coordinate system to a complex (zi,r) coordinate system. 
The inverse, differential and inverse differential forms of the transformation are 
also presented. 

The (zi,r) coordinate system contains three complex planes z\, zi and r. 
The Z3 coordinate is pure real. In order for complex derivatives of a continuously 
differcntiable function f(zi,r) to be holomorphic, the value of the derivative 
cannot depend on the direction of differentiation in each complex plane. The 
set of Cauchy-Riemann equations expressing this condition is determined in 
section 3. 

One of the interesting features of the (zi,r) coordinate system is that the 
XiX2-plane is decomposed into two complex planes z\ and z%. This means that 
each [x\,X2) point in the real plane has one image point in the z\ plane and 
a second in the Z2 plane. The reason for this is that owing to the constraint 
of the Cauchy-Riemann equations not every function f(x\,X2) can be usefully 
represented in a single complex plane. It will be shown in section 3 that the 
transform f(z\,Z2) of any function f{x\,X2) will be holomorphic in both the 
Z!-plane with z% held constant and the ^-plane with Z\ held constant. This 
further implies that any function of N real coordinates will have an holomorphic 
representation in complex space consisting of N complex planes providing both 
spaces are related through an isometric transformation. 

The hydrogen atom is of interest in this paper for two reasons. One is 
that Schrodinger equation generates closed form analytical solutions for the 
hydrogen atom in spherical polar coordinates. The other is that the hydrogenic 
wavefunctions can have angular momentum. By contrast, it is known that a 
field over a one complex plane must be irrotational. It is therefore a test of the 
aforementioned decomposition principle to see if the problem of the hydrogen 
atom can be represented across multiple complex planes. 

In section 4, the Schrodinger equation of the hydrogen atom is transformed 
into both complex Cartesian and complex spherical polar coordinates. One 
advantage of complex versus real spherical polar coordinates is that there are 
no trigonometric functions in them leading to a simpler Laplacian term in the 
Schrodinger equation. It is shown the use of complex spherical polar coordinates 
to eliminate trigonometric functions alongside an application of the complex 
time coordinate to remove the exponential term in the radial component of the 
hydrogenic wavefunctions leads to a much simplified solution process in complex 
space compared to real space. It is further shown that any of the results obtained 
in complex space are readily transformed back into real space. In particular, 
the hydrogenic wavefunctions obtained in complex space are found to transform 
into the generally accepted results [3] in real spherical polar coordinates. 
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2 Complex Space 



The task ahead is to present an isometric conformal transformation relating 
a real (xi,t) coordinate system (i — 1,2,3) and a complex (zi,r) coordinate 
system. This mapping is to be applied in section 4 to a hydrogen atom consisting 
of an electron of total energy E at a separation r(= \xi\) from the proton. It 
can be expressed in the form 

Xi+lX 2 X 2 +IX\ ... 

Z3=X 3 , T = t-l— (2) 

where a is a constant and fi is Planck's constant divided by 2ir. It will further be 
assumed that the components of the , r) coordinate system are independent 
such that 

dzj dzi 

where dij is the Kronecker delta. It is important to be clear that every point in 
the (xi, X2) plane has an image in both the z\ and z-i complex planes. Eq. (JTJ 
therefore constitutes a decomposition of a real plane into two complex planes. 

In the (xi , ^-coordinate system, the electron has a spatial displacement Xi 
from the proton but shares the same world time t. Similarly, in the (zi,r)- 
coordinate system, the electron has a displacement Zi = Xi from the proton but 
shares the same complex time r. It can be shown using eqs. (JlJ and ^ that 

z\z x + z%z 2 + z^z 3 = x\ + x\ + x\ (4) 

(5) 



£1 _ £2 _ 



— 2iz\Zi = xi + x\ (6) 

where z* denotes the complex conjugate of Eq. Q alongside the under- 
standing that the electron and proton share the same value of t confirms the 
isometric nature of the transformation. 

The inverses of the transformation equations (JTJ) can be written 

zi - iz 2 z 2 - izi 

ft \Zi\ 

Xi = Zi, t = r + i— (8) 

hj a 

where \zi\ = \x t \ 

Derivatives with respect to (zi,r) coordinates can be evaluated using the 
chain rule for partial differentiation to give 



d _ 1 
~dz~i ~ 



d d (x\ — 1x2) ifi d 



dx\ dx 2 a\xi\ E dt 



0) 
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d J_ 

dz 2 ~ V2 
d _ 
dz 3 



d 
dx 2 

d 

8x3 ' a\xi \ E dt 



d (x 2 — ix\) ih d 



dxi 

X3 ih d 



a\xi\ 

d_ 

dr 



E dt 

d_ 

~~ dt 



(10) 
(11) 



having used eqs. ([3]) and ([7]). Eq. (JT|) can be used alongside these results to 
give 



_d_ 

''dzi 



d 

' dxi 



\xi \ ih d 
~~E~di 



(12) 



These result will prove uselul later. 
The inverses of eqs. (HJ and (fTO)) are 
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(z\ — iz 2 ) 1T1 d 


dx\ 
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dz 2 


a\zi\ E dr 
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(z 2 — iz\) 1T1 d 


dx 2 




_dz 2 


dzi 


a\zi\ E dr 
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Z3 ih 


d d d 



8x3 



dz 3 



\z,\Edr dt dr 



(13) 

(14) 
(15) 



Eqs. dU) through (fTTj) can be used to validate eq. (]3j) confirming the (zi,r) 
coordinates are all independent of each other. They also give 



dz l 
dz* 



dr 

dr* 







(16) 



showing complex conjugate coordinates are independent of each other. 

It is important to be mindful of the difference between the complex conjugate 
of a derivative and derivatives with respect to complex conjugate variables. For 
example, the complex conjugate of eq. © is 



d* 
dz\ 



1 
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d 
dx\ 



d (x\ + ix 2 ) ih d 



(17) 



dx 2 a\xi\ E dt 
By contrast, it can be shown using the chain rule of partial differentiation that 

(xi + ix 2 ) ih d 



d_ 

Ih* 



1 

V2 



d d 
dxi dx 2 



a\xi\ E dt 



(18) 



Note the sign difference between these two results. The conjugate derivatives 
used in this paper will all be of the form d* jdzi. Conjugate derivatives of the 
form d/dz* have been used in previous work both to represent ladder operators 
in the treatment of harmonic oscillators [5] and to demonstrate the formal elim- 
ination of potential terms from quantum mechanical equations [6] (compare to 

0)- 
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3 The Cauchy-Riemann Equations 



The purpose of this section is to follow through the consequences of the Cauchy- 
Riemann equations for a continuously differentiable function f{xi, t) transformed 
into each of the complex planes z\, z-i and r. These conditions are can be ex- 
pressed as 

df 



df 

at 




dx\ 


ox 2 


dl 

dx 2 


+ i dl 
dx\ 



dr* 

dl 

dz{ 

dl 

dz* 



= 0, (19) 
= 0, (20) 

T, Z 2 ,Z 3 

= 0, (21) 

T,Z1,Z3 



having set 



(22) 



Tir 

S ~ ~oE 

such that t = t + is. Here, the complex variables outside the brackets must be 
held constant during the differentiation. 
Eqs. (jT9|) can also be written 

2 f a 2 E 2 d 2 f \ 

^ + —^ = Vzzz = ^ (23) 

' zi,z 2 ,z s 

It simplifies slightly if / can be assumed to take the product form 

f(Zi,T) = f z (Zi)f T (T) (24) 

to give 

though eq. ([23]) can also be evaluated directly and a product form is not a 
requirement. Any function / will satisfy eq. (|20l) in the form 



dl 

dz\ 



(26) 



since z\ and z* are independent and all of z% , z% and r are held constant. 
Similarly, any function / will satisfy eq. (|2ip in the form 



dl 

dz* 2 



= 0, (27) 

Zl,Z 3 ,T 



since zi and z\ are independent and all of z\, z% and r are held constant. 
It is concluded on this basis that any continuously differentiable function / 
transformed from (xi, £)-space to (zi, r)-space will be holomorphic in the (zi, r)- 
space providing it satisfies eq. 
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Derivatives with respect to real (xi, ^-coordinates that are constrained to 
complex planes can be evaluated using eqs. (fl"3]l through (fT5j) to give 



/ nm+n f \ / i a \ m+n 

( J I / , L (2g) 



V dxfdx% ) ZuZ3tT ~" \.V2dz 2 

df\ df fdf\ df 



— — ) /■ ( 29 ) 

(30) 



dx 3 J , , _ 9^3 ' \dt ) , , , dr 

' 2i,22,T v ' ^1,2:2,23 

where m, n = 0,1,2.... Eqs. ([28)) and (|29]l show that any function f z (%i) will 
satisfy Laplace's equation 



®A ' : 



(31) 



in both the Z\ and z 2 complex planes irrespective of the form of /. This is, in 
fact, just an alternative statement of the Cauchy-Riemann equations (|26[) and 

em. 



4 The Hydrogen Atom 

The Schrodinger equation determining the wavefunction ^{xi,t) for an electron 
of mass m e bound in the Coulomb field of a proton can be expressed in the form 



2m e 



-V 2 * 



47renr 



(32) 



In this, V 2 (= <9 2 /9x 2 ) is the Laplacian operator, -e is the charge on the electron, 
eg is the permittivity of free space and 



E^ = ih 



<9* 

at 



(33) 



gives the total energy of the electron. 

In developing the connection between complex (z iy r)-coordinates and the 
hydrogen atom, eqs. © through (fTTj) and (|33| can be combined to give 



d 
d 

8Z2 



1 
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V2 



d 

dxi 



d (xi — 1x2) 



3X2 

d 



a\xi\ 
(x 2 - IXi) 



a 

dx 2 8x2 ' a|xi| 

5 d £3 

a^3 9x3 alxil 



(34) 

(35) 
(36) 
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These results lead to the operator relationship 



d d* f 8\ 2 d 2 1 

2 = "a — I a — ) z * « h ~2 ( 37 ) 

aa;^ azi \ az^ / ar ozi ar a 

having also used eq. (|T2"j) . It will be convenient to define the V 2 operator in 
complex space to be 

This definition recognises the additional complex terms in eq. (|12p are uniquely 
attributable to the use of a complex time and would vanish if r were reset to t. 
Eqs. (fTTj) and (f37| enable eqs. (|32|) and (f33|) to be rewritten as 

ft 2 T ft 9* 1 / e 2 ft 2 



V 2 * + «(- Ur = (39) 



2m e m e ar dzi r \ Ane m e a 

#7 



E^> = i— (40) 



having set 



2m P E 



(41) 



Eqs. ([55)1 and (HTj) together constitute a description of the hydrogen atom in 
terms of (z{, r)-coordinates equivalent to eqs. (|3"2")l and (f3"3")l . It is usual to solve 
eq. (|32|) using real spherical polar coordinates {r,9,<j>). Eq. (f39|) could also be 
solved in this coordinate system. In this case, we would set 

— - — 

1 dzi dr 

and write V 2 in standard (r, ^-coordinates. The resulting solution would be 
equivalent to the solution of eq. (|32[) but expressed in terms of the complex 
time coordinate r instead of the real one t. 

For the solution of the Schrodinger equation (139]) in this paper it is intended 
to use complex spherical coordinates (z r , zg, z^). These have been chosen to take 
advantage of spherical symmetry but avoid the use of trigonometric functions. 
These complex spherical polar coordinates are defined to be 

Z r =\Zi\, Z0 = -r-^r, z$ = — ^ (43) 

They are related to spherical coordinates in real space through the following 
expressions: 

r = z r , 9 = arccos(ze), 4> — — i\n(Ztj,) (44) 

Eqs. (l44l) are easily checked using eqs. ^ and ([2]) alongside the relationship 
z\ = \zi\exp(i<f>). 
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The calculation of the V 2 -operator in complex spherical polar coordinates 
is similar to its calculation in real spherical polar coordinates but simpler. It 
takes the form 



y2 = -L 



,_d_ 

dz r 



L 2 

h 2 z 2 



where 



L 



(i -4) 



(i 



dz e 



(45) 



(46) 



is read off as the quantum mechanical orbital angular momentum squared oper- 
ator. One of the easiest checks on cq. (j4"5|) is to transform it in to real spherical 
polar coordinates using eqs. (|4"4"|) . 

It will be convenient to assume the Schrodinger equation (j3"9")l has a product 
solution of the form 



*(z r ,z e ,Z0,r) = R(z r )@(z e )G>(z<f,)exp(-iET/h) 
enabling it to be decomposed in to the component equations 



d 2 R 



dR 



+ ( 2 " 0-77 — I 1 — A + -j- ) R = 



d£ 

d 

dze 



dze 



L 



fc 2 9 



L 2 Q 



= 



where 



Zcj, 

2r 
a 



dz d 



$- fc 2 $ = 



A = 



—m 
~2E 



and L and k are the constants of separation. 

It is notable that L and k must satisfy the eigenvalue equations: 

= LH, 

where the form of the z-component of angular momentum operator 

f 9 

OZa, 



(47) 

(48) 
(49) 
(50) 

(51) 

(52) 
(53) 

(54) 



can be identified from eq. (|5D|) . 

Eq. (|5"3")) has the eigensolution $ = z\. This transforms back to $ = 
exp(—tk(j)) in real spherical polar coordinates. Clearly, k must be and integer 
for $ to be a periodic function. 

Eq. (|49|) is Legendre's equation [8] giving 6 = Pj°(zo) where P? are asso- 
ciated Legendre polynomials, L 2 = 1(1 — \)Ti 2 and Z(> 0) is the orbital angular 
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momentum quantum number. In this case, the connection to real spherical polar 
coordinates is immediate since zg = cos(#). 

The solution of eq. ([18} is facilitated through the transformation R = t; l R' 
into Laguerre's differential equation. Putting these results together leads to the 
complete solution of the Schrodinger equation for the hydrogen atom in complex 
(z r , zg, z<p, r)-coordinates. This takes the form 

=Kl (—) l L^+l_ 1 (^) Pfize^eM^Er/h) (55) 
\na J \na J v 

where M n i is the normalizing constant and the integer n{> 0) is the principal 
quantum number. The energy eigenvalues for the bound electron 

4 

„ m e e 

E n = 56 

are determined during the solution of the radial component of the equation. 

Eq. (|55[) is readily transformed into (r, 9, <fi, r)-coordinates to reveal the form 
of the hydrogenic wavefunctions most often presented in standard textbooks on 
quantum mechanics. It is notable, however, that the solution of eqs. (|48p 
through (|50p has taken fewer steps than the solution starting from the initial 
component form of eq. (|32p in (r, 9, <fi, t)-coordinates. This is for two reasons. 
First, the complex time coordinate served a useful purpose in eliminating an 
exponential term of the form exp(— 2r/an) in the radial wavefunction. This is 
otherwise a manual process in (r, 9, <f>, t)-coordinates in order to reduce the radial 
problem to Laguerres's equation. Second, the use of complex versus real spher- 
ical polar coordinates has stripped away the need to manipulate trigonometric 
functions. 



5 Concluding Remarks 

The motivation for this paper has been the idea that problems that cannot 
be represented in one complex plane can be decomposed for representation over 
multiple complex planes. There are two reasons a problem might require decom- 
position. One is that it generates fields that do not satisfy the Cauchy-Riemann 
equations. The other is that the problem is defined in more than two dimensions. 

The method that has been found for decomposing problems is to use an 
isometric transformation to take the problem from a real vector space to a 
complex one. This term decomposition applies since any point (2:1,2:2) in a 
plane in the real vector space will transform to (z\,z-i) and therefore image 
points in two complex planes. It has been further shown any function j(x\,X2) 
will be holomorphic in both the zi-plane with zi held constant and the Z2-plane 
with z\ held constant. It is inferred that any problem that can be expressed 
in a real vector space of iV-dimensions can also be represented in a complex 
space composed of no more than N complex planes such that all fields will be 
holomorphic in the complex planes. 
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The decomposition principle just outlined has been demonstrated for the 
case of the Schrodinger equation of the hydrogen atom. This well known sys- 
tem has both more than two dimensions and features such as angular momentum 
that cannot be represented in a single complex plane. It has been nevertheless 
transformed into a complex space consisting of three complex planes and solved 
for the complete set of hydrogenic wavefunctions. The solution has utilized com- 
plex spherical coordinates that eliminate the need for trigonometric functions. 
The Laplacian operator has been found to be simpler in complex versus real 
spherical coordinates facilitating a solution process that overall has many fewer 
steps in complex space than in real spherical polar coordinates. 
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